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NONLINEAR INITIALIZATION OF THE GLAS MODEL 
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ABSTRACT 

A simplified version of the GLAS model is linearized 
and the normal modes extracted. These modes show the 
necessary separation for nonlinear initialization. 

INTRODUCTION 

The finite difference form of the GLAS model is exceptionally 
difficult to put into normal mode form. A procedure is outlined 
which will allow a simplification of the model and still incorpo- 
rate the essential ingredients for nonlinear initialization so 
that the high frequency components may be slowed. A small version 
of the simplified system is presented, linearized, and the normal 
modes described. Some form of this system will ultimately be 
initialized and tested in the GLAS model. 

INITIALIZATION 

The procedure for nonlinear initialization has been presented 
by Baer and Tribbia (1977). The system equations should be 
written in the form. 


+ A X = e G(x,x) (1) 

where X is a vector of dependent variables (grid point values of 
the flow field and physical variables, for example), A is a matrix 
of fixed quantities dependent on the system, and e issmall paran- 
eter, usually the Rossby number. Since most systems of this sort 
are exceptionally complicated, as in the GLAS model, let us re- 
write (1) in the form 


+ A X = e G + (A- A) X (2) 

where A is a "simpler" matrix, and || S _1 (A-A)s|| is of order e. 

Here S is the modal matrix which^diagonalyses A. 

We shall find the modes of A, a simple version of the GLAS 
model so that we may separate the frequencies into fast and slow, 
a procedure essential to the nonlinear initialization scheme. Thus 
we will find, 
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where clearly the scaled frequencies eX << X 

x y 

GLAS MODEL NORMAL MODES 

We begin the simplification of the GLAS model from the non- 
forced version (see Somerville, et. al. , 1974) described as 
follows : 
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Symbol definitions may be found in the above reference. We now 
linearize this system about a state of rest, let n -* Tr(x,y,t) + II 
and a -* a(x,y,a,t) + A(a) where tt << IT and a << A. We include the 
following definition where <J> is the perturbation of $ , 

* = * + oAir,u = o+ £|J,V = (u,v) . (5) 

The resulting linearized system then becomes, 
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If we now assume a solution in longitude of periodic form, 
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Now to separate the vertical dependence from the latitudinal, 
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(u,v,t|0 = G(o) (u,v,^) (y , t) 
u = H(a) u(y,t) 

Substitution of (8) into (7) yields the result that 
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and application of the boundary conditions that o = 0 at a = 0,1 , 
gives the conditions that H(a=0) =* 0 and H (a=l) = - (TI | C^)A(l) 

H (1) . These conditions applied to the solution of (10) yield 
the required separation constants C 2 , which are also denoted as 
"equivalent depths". 

The first two equations of (7), together with (9) allow for 
a solution of the latitudinal structure and generate the normal 
modes and frequencies. We have, 
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To better understand these modes in a simplified framework, let us 
solve (11) in a channel extending from - 60° £ 6 £ 60° latitude 
with rigid walls so that v(+60) = 0. Equations (11) also show 
that v yy (+60) = 0 and give equations on the boundaries for and 
u. For additional simplification, we assume the channel to be 
broken into 2N + 1 equally spaced grid-points such that Y = NAy 
and - Y £ y < Y, with y = nAy. We also use centered differences 
to replace derivatives. In the finite-difference form, system 
(11) becomes, 

' |*+4x = 0 x T = (U.V,*) 
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and A is a 3(2N-1) square matrix. The roots of A represent the 
frequencies of the reduced system as specified in (3) and will 
represent an approximation of the GLAS model. 


STRUCTURE OF MODES AND FREQUENCIES 

For pilot experiments, we have chosen to break up the atmos- 
phere into three equal a-layers. Finite difference solution of H 
from (10) then yields the following three equivalent depths in 
units of geopotential height; 
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= - 8.402*10 4 
C 2 2 = - 4.235*10 3 
C 2 = - 8 . 012 * 10 2 

The corresponding structures of vectors H ,H and may be seen 
on Figure 1. They have the expected distribution for this re- 
solution. For the latitudinal structure, we have taken Ay = 5° 
and Ay = 10°. For the former case, we generate 69 modes for each 
equivalent depth (C^) and for each planetary wave (k) . For the 
latter case, the number reduces to 33 modes. In either case, 
many of the modes are computational and presenting their struct 
tures would be overwhelming. Fig. 2 gives some indication of the 
structure of the latitudinal modes. We have plotted the 49th 
vector for the case of Ay = 5°, showing the distribution for 
waves 1, 6, 15 and 30 for both the external and first internal 
modes. Note that the vector for wave one and the external mode 
shows some truncation properties which are not apparent for the 
internal mode. The other vectors all have a very similar struc- 
ture. These vectors reflect Rossby type oscillations with 
periods in excess of 40 hours. 

In Table 1 we describe the distribution of eigenvalues for 
the Ay = 10° truncation. We denote the number of modes which fall 
into the n high M frequency as compared to the "low" frequency for 
selected wave numbers and the three vertical structures. Note 
that all but wave number one have very strong separation between 
the gravity and rotational modes; i.e., there are twenty two "fast 1 
modes and eleven "slow" ones. Separation in wave one, although 
not pronounced, is apparent. This difficulty with long waves has 
been seen in other experiments. 

SUMMARY 

The data presented indicate that normal modes may be ex- 
tracted from a simplified version of the GLAS model, and mode 
separation is evident to allow for nonlinear initialization. Tests 
must now be performed which will determine how important the 
approximations to the true numerical model are, and how effective 
initialization with the simplified model will be. 

An interesting observation of the study to date shows that 
some of the modes are strongly affected by truncation. Devices 
for removing the amplitudes of these modes (if they have slow 
frequency) are under investigation. The high frequency computa- 
tional modes may be removed by the nonlinear initialization 
technique. 
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TABLE 1: Periods in hours for all frequencies for model with 

Ay = 10° and for selected wave numbers. Values in 
parentheses indicate the number of modes in the 
specified range. 


WAVE NUMBER EXT- MODE FIRST INT . SECOND INT . 


1 

6 < T < 16 (19) 
23 1 T 1 45 (7) 

T > 189 (7) 

15 < T < 56 (22) 
T > 170 (11) 

15 < T < 60 (20) 
T - 102 (2) 
T > 393 (11) 

6 

4.1 < T < 6.5 (22) 
T - 18.2 (2) 
T > 85 (9) 

13 < T < 29 (22) 
T > 103 (11) 

15 < T < 43 (20) 
T « 65.6 (2) 
T > 140 (11) 

15 

2.3 < T < 2.6 (22) 
T - 31-9 (2) 
T > 198 (9) 

8.9 < T < 11.4 (22) 
T > 201 (11) 

13 < T < 26.3 (22) 
T > 238 (11) 

30 

1.24 < T i 1.28 (22) 
T - 58.8 (2) 

5.3 < T < 5.7 (22) 

T > 312 (11) 

10 < T < 13.2 (22) 
T > 420 (11) 
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Figure 1: Vertical structures for the external and two internal modes of the 

simplified model. The vectors are not normalized. 
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